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We study the effects of bounded confidence thresholds and of interaction and external noise on 
Axelrod's model of social influence. Our study is based on a combination of numerical simulations 
and an integration of the mean-field Master equation describing the system in the thermodynamic 
limit. We find that interaction thresholds affect the system only quantitatively, but that they do 
■ not alter the basic phase structure. The known crossover between an ordered and a disordered 

' state in finite systems subject to external noise persists in models with general confidence threshold. 

, Interaction noise here facilitates the dynamics and reduces relaxation times. We also study Axelrod 

. systems with metric features, and point out similarities and differences compared to models with 

nominal features. Metric features are used to demonstrate that a small group of extremists can have 
a significant impact on the opinion dynamics of a population of Axelrod agents. 
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PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



o 
o 

Given the increasing influence of mass media, globalisation, electronic communication and intercontinental travel 
the apparent persistence of cultural diversity seems surprising. To study this problem Axelrod [l[ proposed a simple 
agent-based model to study how cultural features disseminate. In particular the model addresses the question of how 



cultural diversity can result from locally attractive interaction. Axelrod's model is -in the language of physics and 



dynamical systems theory- a cellular automaton with a set of discrete degrees of freedom placed on a discrete spatial 
lattice, updated in time through specific interaction rules, and is easily simulated computationally. The simulation and 
theoretical analysis of social and economic systems has in the recent years been taken up by the statistical physics 
community (see e.g. 0, @, 0, S @ and references therein). Social systems are here interpreted as many-particle 
problems and tools originally developed to study physical systems have been transferred and adapted to the study of 
socio-economic models in different contexts. See also [7| for a sociological perspective of agent-based modelling, 
(■sq One of the first systematic studies of the Axelrod model is the one of [8(, where a phase between a monocultural 

state (referred to as globalisation) and a multi-cultural phase (referred to as polarisation) has been identified as a 
function of the degree of variation of the initial conditions chosen before running the Axelrod automaton. In terms 
\ of physics this transition is one between ordered and disordered states. Axelrod's two-body interaction is attractive 
as such, and drives the system to order. Two particles can only interact, however, if their states are not fully distinct 
[~^. \ from each other. These interaction barriers lead to potential jamming. The model thus exhibits some similarity with 
kinetically constrained models of glass forming materials 0, [l(| ■ If the disorder of the initial conditions is higher than 
some critical value, the kinetic arrest occurs after a finite number of time-steps, and the systems remains jammed in 
a disordered state. For low enough initial disorder in the random initial conditions, however, the ordering dynamics 
continues indefinitely in infinite systems. In finite systems full convergence is reached at finite times, and the dynamics 
comes to a halt when any remaining disorder has been eliminated, resulting in a fully monocultural 'globalised' state. 

A variety of extensions and variations of Axelrod's original model have been proposed, Klemm et al. have for 
example studied the effects of noise on the jamming behaviour of Axelrod systems [111 . Il2j |. furthermore the effects 
of mass media influence has recently been addressed in [l3|, The model, originally defined on a square lattice, 
has furthermore been simulated on a variety of complex networks in order to study the effects of the topology of the 
underlying web of interactions on the dynamics and convergence properties [l5j . 

One of the obvious shortcomings of Axelrod's original model and the subsequent variations is the fact that a metric 
structure is missing in opinion space. Opinions Oif of agent i on feature / take values in a discrete set, usually labelled 
by {1, . . . . , q}, and agents are classified only as to whether they have the same opinion on a certain issue (<7j/ = Oj/), 
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or whether their opinions are different (aif ^ a jf)- No notion of partial agreement on a feature is present, agents can 
either fully agree on a certain issue / or they disagree (a more formal definition of the model will follow below). This 
assumption has been relaxed in a sociological context for example in [lfl [13, EH , and the notion of metric spaces has 
been introduced. For such features a gradual distinction of a degree of agreement can be defined, for example given 
by |er,;/ — Ojf\ m , where the integer m is referred to as 'moderation' (l6| . 

As a second drawback, no notion of different tolerance levels or inclination to change one's opinion is present in the 
original setup. Interaction between agents in Axelrod's original model can occur once they agree on at least one out of 
a number F of features. F along with q mentioned above are the main model parameters in Axelrod's original setup. 
From a sociological point of view it is interesting to study the model in more generality, and to introduce a 'confidence 
threshold' # € {0, . . . , F}, so that agents have the potential to interact only if they agree on (strictly) more than 1? 
out of F issues. This is referred to as 'bounded confidence' in the sociological literature, Axelrod's original model 
corresponds to minimal confidence threshold (i? — 0), in which interaction is possible rather easily. Large values of the 
threshold d systematically suppress the potential to interact, corresponding to more and more conservative agents, 
who do not change opinion easily. 

Studies of different models addressin g eit her of these two points e.g. through the introduction of continuous opinion 
states can be found in [ijj [2(| [U I22TI231 [24], [25| . A first analysis of the effects of metric features and confidence 
thresholds in the context of the Axelrod model has been conducted in [13, OH- These studies focus mostly on numerical 
simulations and is mostly restricted to a specific choices of the model parameters F and q. The aim of the current 
work is to complement and extend the analysis of [171 . 1 1 81 ] through a more general study of the model in parameter 
space. We also provide analytical results based on a Master equation approach 0] for the Axelrod model with nominal 
features and general confidence threshold. In addition, the introduction of a metric model allows us to address issues 
such as extremism in the context of Axelrod opinion dynamics, which are not captured by the conventional 'nominal' 
formulation. 



II. THE MODEL 

The system is composed of N = L x L agents fixed on the nodes of square lattice of lateral extension L. For 
simplicity we consider periodic boundary conditions in both spatial dimensions. The state of agent i 6 {1, . . . , N} 
is characterised by an opinion vector Bi = (era, ■ ■ ■ , c^f), where the integer F > 1 denotes the number of cultural 
'features' in the model. Each component aif then indicates the opinion of agent i on issue /. In Axelrod's original 
formulation each component aif takes one of the q values {1, . . . , q} at each time step, so that each spin Bi describes 
one of q F cultures. Initially each aif is drawn at random from the set {1, . . . , q} with no correlations between agents 
or features. The model parameter q hence measures the degree of disorder in the random initial spin configuration. 

We will in the following distinguish between 'nominal' and 'metric' features as suggested in [13, Oil]. We first 
describe the dynamics of the nominal Axelrod model. Here, the system evolves in time by iteration of the following 
steps: 

(i) Select one spin i S {1, . . . , N} at random. Subsequently select one of its four nearest neighbours at random. 
Call this second spin j. 

(ii) Compute the overlap uj(i,j) = X)/=i ^o-./^i/ ^ {0, ■ ■ ■ , F} between spins i and j (with S a ^ a i the Kronecker 
delta). 

(iii) If uj(i, j) = F continue with (v) (spins Bi and Bj are in identical states). 

(iv) Set the probability for i and j to interact to I = 8 if uj(i,j) < d and to / = uj(i,j)/F if j) > d. Then 
with probability I — I leave spins i and j unchanged. With probability / spins i and j perform the following 
interaction: choose one feature / at random such that a t f ^ ajf. Such a feature exists as uj(i, j) < F. Then set 

°i/ = a :, '■ ■ 

(v) External noise. With probability 7 perform the following: choose one spin i and one feature / at random. Set 
aif to a value chosen randomly from {1, . . . , q}. 

(vi) Resume at (i). 

We will refer to one cycle (i)-(vi) as a microscopic time-step in the following. At system size N the duration of 
such a step is taken to be At = l/N. In general we will present the time-evolution of the system mostly in terms of 
macroscopic time units i, so that one unit of time t corresponds to N microscopic interaction cycles, i.e. on average 
to one (attempted) update per spin. 
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In the above dynamics $ is the interaction threshold mentioned above. In the absence of interaction noise (S = 0) 
neighbouring agents have the potential to interact if and only if they share opinions on (strictly) more than ■& out 
of F features. To soften this constraint we follow [l7| and introduce a source of noise, and allow agents who agree 
on or fewer features to interact with probability 6. We refer to this type of stochasticity as 'interaction noise' in 
the following, 5 measures its strength. 7 in the above update rules instead denotes the strength of what we will call 
'external noise'. After each time step, with probability 7 a randomly chosen component of a randomly chosen spin is 
set to a random value {!, ... ,q}. This type of noise has first been studied in [Tlj . 



III. MASTER EQUATION IN THE MEAN FIELD APPROXIMATION 

In this section we will consider a mean field approximation of the model. In the mean field model it is possible and 
convenient to consider the dynamics in terms of bonds, i.e. of pairs of neighbouring agents, rather than in terms of 
spins {<Ti}. Following the strategy of [1, US HH , let P m {t) be the probability that, at a given time t, a bond is of type 
m, i.e. that the two agents at the ends of the bond have the same opinion on exactly m features. We will occasionally 
refer to bonds of type F as 'fully saturated' in the following. If we let p be the probability that at the starting point 
of the dynamics two spins share a given feature, we have initially 



F-m 



(1) 



For (Tif drawn independently and with equal probabilities from {1, . . . , q} one has p = l/q. 

We further define A to be the probability that two independent spin components are equal but different from a given 
third. A is in principle a time-dependent quantity as the system evolves according the Axelrod dynamics. We here 
neglect this time-dependence and assume that A is well approximated by its initial value A = (q — throughout the 

dynamics. This was seen not to have any significant effects on results in UHjJ- Now denote by W^m the transition 
probability that a bond of ty pe n becomes of type m due to the updating of a neighboring bond of type k. The only 



non-zero elements are [g, [27], [28 1 



- n/F, 
Wi% = (1 - A) (1 - n/F) , 
= X(l-n/F) , 

independently of k. We will therefore suppress the superscript k in the following. 

Let us further define Ik to be the probability with which two agents who share opinions on k features interact if 
selected for potential update. Then one has 

k<$ 
d+l<k<F 
k = F. 

The master equation can then be written in the form 
gdP m (t) ^! , 




2 dt 



fe=0 
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(1 - S mJ )P m+1 r -^±l (l - i) - (1 - 6 m<f )P m ^- q 
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(2) 



This equation is an approximation in the mean-field sense, and the thermodynamic limit is implied. The Master 
equation can be expected to describe the system at most at large system sizes, and will therefore not be able to 
capture features characteristic of finite systems. The geometry of the square lattice is mimicked, in the mean-field 
spirit, by the pre-factors g — 1 and g in the different terms of the Master equation, g here denotes the co-ordination 
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FIG. 1: (Colour on-line) Relative size S/N of the largest cultural region in the standard Axelrod model with F nominal 
features, initial conditions drawn from (Tj/ G {1, ...,<?} with flat distribution. Symbols show data from simulations run until 
convergence, system size N = 400, averaged over 10 samples, F = 5, 10, 20. The vertical dashed lines mark the location of the 
phase transition as predicted by a numerical integration of the Master equation. The right panel shows Pf from simulations 
(symbols) compared to results predicted by the Master equation (solid lines). 



number of each spin so that these coefficients reflect the number of spins with whom a given spin can interact (albeit 
these are not nearest neighbors any longer). On a square lattice in two dimensions one has g = 4. The pre- factor 
g/2 in front of the time derivative in Eq. © takes into account the fact that the system contains g/2 bonds per 
lattice site. One would expect the Master equation to be accurate in the case of degree-regular graphs (of connectivity 
g), as discussed for example in [28(. Still, as demonstrated [8( and as we will see below in the context of external 
and interaction noise, an approach based on numerical integration of the Master equation is able to reproduce some 
features of the two-dimensional model at least qualitatively. 



IV. AXELROD DYNAMICS WITH NOMINAL FEATURES 



A. Baseline model 



For completeness we re-itcratc the behaviour of the baseline Axelrod model (9 = 7 = 8 = 0) in Fig. \T\ For any 
given number F > 2 of features a discontinuous transition between an ordered state at q < q c (F) and a disordered 
phase at larger values of q is observed. At q < q c the coarsening dynamics of the model persists until a fully ordered 
state is reached. For any feature / = 1, . . . , F all agents then agree on one opinion, i.e. <Hf = <jjf for all i,j. In 
finite systems such a state is reached after a finite time. Fig. Q] depicts the relative size S/N of the largest culturally 
homogeneous region of spins as a function of q. A homogeneous region R is here defined as a subset of the L x L 
agents, so that within R all agents agree on all features [29]. As seen in the figure, one finds only one region at 
convergence for q < q c , and has S/N = 1. At values of q larger than q c (F) dynamic arrest occurs before the system 
can reach a fully ordered phase. After the arrest no further ordering is possible due to the kinetic constraints imposed 
on the otherwise attractive spin-dynamics. The system remains in a disordered state, marked by a large number of 
small cultural regions and a vanishing number of active bonds. S/N remains small at convergence in this regime. As 
seen in Fig. Q] the fraction Pp of fully saturated bonds behaves similarly to S/N at convergence, and can be well 
captured by the Master equation in the disordered regime. The ordering at low values of q can not be obtained from 
an approach based on the Master equation [f| . 



B. Effects of interaction threshold 



We now turn to the noise-free model with general interaction threshold As shown in Fig. [2] the qualitative 
behaviour of the model is not affected by the introduction of a confidence threshold. As before an ordered phase 
is found at low values of q, and a disordered one at q larger than some critical value q c ("d,F). One finds that an 
increased confidence threshold suppresses interaction and hence reduces the range of q in which order can be reached. 
q c is a decreasing function of ■& at fixed F. Fig. [2] also demonstrates that the Master equation given above describes 
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FIG. 2: (Colour on-line) Fraction of fully saturated bonds Pf in the Axelrod model with F = 10 nominal features and with 
confidence threshold i? = 3,2,1,0 (from left to right), initial conditions drawn from a%f G {l,...,q} with flat distribution. 
Connected markers are from simulations with N = 400 agents, run until convergence, averaged over 10 samples. Thick solid 
lines show theoretical predictions by the Master equation. 



the qualitative behaviour of the system and dependence on the interaction threshold appropriately. In the disordered 
phase even a reasonable quantitative agreement between numerical measurements of the fraction of fully saturated 
bonds and the corresponding theoretical predictions can be observed. We attribute remaining discrepancies to the 
mean field approximation, inaccuracies in capturing the 2-dimensional geometry and to finite-size effects. 

Fig. [3] depicts the phase diagram of the model in the q) plane for different values of F, as obtained from the 
Master equation (jjjj. The disordered phase is found at large values of $ and q respectively, order is reached at low 
q and/or It may here be interesting to ask whether the relevant variable is the absolute interaction threshold 
$ S {0, . . . , F}, or the relative one d/F. In [l?], [H| results are for example reported in terms of relative thresholds. 
The left inset of Fig. [3] confirms that the phase boundaries for different values of F as shown in the main panel do 
indeed show a reasonable collapse if plotted as a function of •d/F, At small values of $ systematic deviations are 
however observed. A different rescaling was suggested in [3l|, where results for the one-dimensional Axelrod model 
where shown to depend mostly on q/F. As demonstrated in the right inset of Fig. [3] equally good collapse is observed 
in the (i!}/F,q) plane, so that we can here not reach a definitive conclusion as to whether there are any independent 
scaling parameters, and if so which ones they are (32j. 

Some indications regarding the relevance of absolute as opposed to relative thresholds can be found in Fig. 2J 
where we show the density of fully saturated bonds Pp as a function of the density of initially active bonds n a (0) = 
Sfc=i9+i Pk(t — 0). Simulations are here performed by fixing $ and F and by subsequently varying q. n a (0) then 
decreases with increasing q. The data shown in the figure suggest a potential collapse on three different curves, one 
for each of the tested values i? = 0, 2, 8. While these findings might indicate some potential universality as F and q 
are varied at fixed reaching a final conclusion as to whether absolute or relative thresholds are the relevant ones 
still remains an open question. 



C. Effects of external noise 



We now turn to a discussion of the effects of noise on nominal Axelrod dynamics. We will first focus on external 
noise as introduced above. This type of stochasticity describes random fluctuations which are triggered by an external 
event, and which result in individual spin components being flipped randomly at a given rate. External noise was 



introduced in the context of the Axelrod model in 11(, and the resulting random mutations describe what Axelrod 



refers to as 'cultural drift' [l| in the population of agents. Klemm et al. [TT| have studied the effects of external 
noise as a function of the noise rate and system size in fully equilibrated systems. We here extend this analysis, and 
consider different cases according to whether the thermodynamic limit or long-time limit is taken first. Results from 
an integration of the Master equation are discussed in order to provide a semi-analytical description of the system in 
the limit of infinite size. 
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FIG. 3: (Colour on-line) Phase diagram of the model with general confidence threshold in the mean-field approximation. 
Locations of the transitions are obtained from the Master equation. The curves show F — 40, 20, 10, 5 from top to bottom. 
The system is in the disordered phase above the respective lines. Ordered states can be expected below. The insets show a 
rescaling in terms of the relative threshold ■d/F (left) and of the relative number of opinion states per feature q/F (right). 




FIG. 4: (Colour on-line) Pf versus fraction of active links n a (0) in initial configuration for nominal Axelrod dynamics at 
different & and F. Data are from simulations (system size N — 20 x 20, run until convergence, averages over 10 samples are 
taken). 



1. Finite system, equilibrated dynamics 



In finite systems a continuous transition between an ordered state at low noise rates 7 and a disordered state at 
large 7 has been identified in This transition relates to a characteristic relaxation time T = 0(N log N) in finite 
systems. When the noise rate is sufficiently large (larger than T _1 ) stochastic perturbations build up in time, and 
lead to a disordered state. For 7T << 1 the system drifts from one ordered state to another in time, time-averaging 
effectively yields global order. Since at q > q c disorder is observed in the absence of noise (7 = 0), the behaviour of 
the model in the limit 7 — > + is discontinuous at q > q c . As discussed in the next point the finiteness of the system 
is crucial here, so that the described behaviour cannot be captured by the Master equation. 

The effects of external noise on finite Axelrod systems with general threshold is depicted in Fig. [5j The behaviour 
of the model with confidence threshold is here found to be very similar to the one identified in for conventional 
Axelrod dynamics. We here focus on F = 10, q — 100 as an example, but similar behaviour can be expected for 
other model parameters in the disordered phase of the noise-free model. For small values of 7 the system orders 
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FIG. 5: (Colour on-line) S/N versus time for F — 10, q = 100. # = 0, 1, 2 in top, middle and lower panel. Curves shown for 
7 = 10~ 2 , 10~ 3 , 10 -4 (bottom to top at large times). From simulations, N = 10 X 10, averages over 10 runs. The a:-axes show 
time in multiples of N microscopic updates, i.e. after t every agent has on average been selected f-times for potential update. 

after an initial transient. At large magnitude of the applied noise, no ordering is found, consistently with the results 
of [Tl| . This general qualitative picture appears to be independent of the applied threshold. The duration of the 
equilibration period in cases where the system orders, however, shows a significant dependence on the noise strength 
and on the chosen threshold. Generally, the time required to reach equilibration increases as 7 is lowered or as $ 
is increased, see Fig. [5j The value of S/N at equilibrium is a decreasing function of 7. In the examples of Fig. [5] 
equilibration to a value of S/N « 0.1 occurs fast at 7 = 10~ 2 . At lower noise rates, S/N reaches values in the range 
of 0.8 to 1, but only after a substantial equilibration period, which increases as 7 is lowered. Only models with low or 
moderate threshold and/or sufficiently large noise strength can hence be equilibrated in reasonable computing time. 
While analogy suggests that an ordered phase sets also at higher thresholds and small enough noise strengths if 
the dynamics is run long enough, we have not been able to confirm this explicitly due to computational limitations. 
Approaches based on continuous-time Monte Carlo methods might here potentially be more appropriate than direct 
simulation of the Axelrod dynamics [33] . 

2. Thermodynamic limit for equilibrated systems 

Due to the divergence of the relaxation time T with N, the transition just discussed disappears in infinite systems. 
If, at q > q c , the system size is taken to infinity after equilibration, i.e. if equilibrated systems are considered at 
increasing N while keeping all other parameters fixed, the population always ends up in a disordered state for all 
7 > 0, as demonstrated in [111 ]. The discontinuity at 7 — > + is removed. Equilibrated systems of infinite size cannot 
be captured by the Master equation, as the latter implies the thermodynamic limit to have been taken first. 

3. Thermodynamic limit at finite t 

Taking the thermodynamic limit N — * 00 at a fixed number of macroscopic time steps t results in non-trivial 
behaviour. As displayed in Fig. [5] one finds a disordered state with Pp w at sufficiently large 7. Partial ordering 



8 




Y 

FIG. 6: (Colour on-line) Pf versus 7 in a nominal Axelrod system, F = 10, q = 100. Upper panel: results from simulations, 
run for t macroscopic time-steps, where t = 20000, 10000, 5000, 1000 from left to right. Circles are for N = 20 x 20, squares 
for N = 30 x 30, diamonds TV = 40 x 40, triangle N = 50 x 50. Data are averages over 10 samples. Lower panel: Pf from 
numerical integration of the Master equation, t = 20000, 10000, 5000, 1000 as above. An Euler-forward scheme with time-step 
dt — 0.1 is here used. Dashed line in both panels marks result from Master equation after even longer times (t — 2 • 10°). 

sets in as 7 is lowered, with non-monotonous behaviour and a peak of Pp at intermediate noise-strengths. As 7 is 
reduced further, Pp takes small but non-zero values, mostly independent of the noise strength, provided the latter is 
small enough. Integration of the Master equation confirms this behaviour qualitatively. The non-monotonic behaviour 
of Pp as a function of 7 at fixed time-scale t is here potentially related to the non-monotonic temporal behaviour of 
the Axelrod system as observed for example in [g, |28| . Under suitable conditions order parameters such as Pp (t) or 
S(t)/N might become non-monotonic functions of time at fixed values of q,F,j. These non-monotonicities are then 
reflected as peaks in Pp when other cuts through parameter space are considered, as in Fig. [5] where {t, F, q\ are 
fixed and 7 is varied, or in Fig. [7] where q is varied at fixed {F, t, 7}. 

4- Large time limit after taking thermodynamic limit 

As the time-scale on which the systems is studied is increased, the peak in Pp appears to move further to the left 
in Fig. [6l i.e. to smaller values of the noise strength 7. Hence for any fixed 70 there is a time-scale £(70) so that P(F) 
is monotonically decreasing as a function of 7 > 70 on this time-scale. The system is hence disordered at large 7, and 
partially ordered at low 7. This behaviour can successfully be described by the Master equation. 

D. Effects of interaction noise 

We next turn to a discussion of the effects of interaction noise, as parametrised by its amplitude 8. This type of 
stochasticity facilitates interaction, as it removes kinetic constraints and allows agents sharing ■& or less opinions to 
interact (at rate 8), while in standard Axelrod dynamics they would not be able to align spin components. While the 
external noise of amplitude 7 has an ambiguous role of inducing order at low amplitudes and of driving the system 
to disorder at large 7, interaction noise can generally be expected to favour order. It his hence interesting to study 
the system in presence of both types of randomness, and to identify re-inforced ordering behaviour or (at small 7) 
potential competition between the ordering and disordering stochasticity (at large 7) 

Results for a nominal Axelrod system with both types of noise are reported in Fig. [5J The data indicate that the 
effects of interaction noise are mostly to facilitate order for large ranges of fixed external noise 7. More specifically, 
the effects of interaction noise is to reduce the time-scale on which finite systems order in the presence of external 
noise. The left panel of Fig. [5] shows the concentration Pp of bonds with full overlap in an Axelrod system run 
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FIG. 7: (Colour on-line) P(F) versus q in a nominal Axelrod system, F = 10. Results are from numerical integration of the 
Master equation up to t = 5000. Curves show 7 = 0.0001, 0.0005, 0.001, 0.005, 0.025, 0.05 from left to right at the maximum. 
Curves for 7 = 0.025 and 7 = 0.05 still display a maximum, but to the right of the plotted range of q. 

for a time which is not long enough for the system without interaction noise (circles, S = 0) to develop order at 
the studied magnitudes 7. Order at 7 < 0.001 would develop only if the system were run for longer times, and had 
fully equilibrated. The curves for non-vanishing amplitude S demonstrate the effect of interaction noise, the system 
now orders at small noise strengths 7 . While the qualitative behaviour of the equilibrated system is not altered, the 
facilitation of the kinetic constraints drastically reduces equilibration times, and the system orders at sufficiently low 
7 even after moderate running times. Interestingly the noise strength 70 separating the ordered from the disordered 
regime of the fully equilibrated system appears not to be affected much by the interaction noise. The curves displayed 
in the left panel of Fig. [8] are indeed mostly independent of S, as long as 5 > 0. It mi ght potentially be interesting to 
study even lower S, although probably unrealistic from the sociological point of view [34j . 




We conclude that the effect of interaction noise is to reduce relaxation times, but that it does not alter the phase 
behaviour of the model, with an ordered phase at low 7, and a disordered one at large 7. As in the absence of 
interaction noise (S = 0) this transition is present only in finite systems, seen in the inset of the left panel of Fig. [8] 
and in the right panel. As the system size is increased at equilibrium the order-disorder crossover moves to smaller 
values of the noise strength 7, and can be expected to be absent in the thermodynamic limit, where only the disordered 
region prevails. Indeed rescaling of the data in the inset of the left panel demonstrates that 'yNlnN is the relevant 
scaling variable, similar to the observations of (llj . At finite running times t the order at low 7 is gradually reduced 
with increased system size and in the thermodynamic limit the system is qualitatively well described by the Master 
equation. 

Fig. [^finally confirms that this behaviour is not limited to the standard Axelrod dynamics with vanishing interaction 
threshold ($ = 0). Interaction noise reduces the time-scale on which the system orders at small 7 also in the model 
with moderate non-zero thresholds, and that the system at d > behaves very much like the one at 1) = 0. 



We now turn to a modification of the Axelrod dynamics in which a metric allowing for gradual notions of agreement 
between agents is introduced. As before, opinions take discrete values a if £ {l,...,q}. However, the 'distance' 
between two spins at and <ij is no longer measured in terms of the number of features on which the corresponding 
agents agree, but based on the following Euclidean distance between the spin- vectors c?j and uf. 




V. 



EUCLIDEAN AXELROD DYNAMICS 




(3) 



Thus the distance between two agents ranges between and 1. It takes the maximal value d = 1 if and only if 
the opinions of the two agents are diametrically opposed, i.e of for any feature / one has a if — l,<Jjf — q or 
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FIG. 8: (Colour on-line) Effects of interaction noise. Left: Pf versus 7 for Axelrod model (F = 10, q = 100, ■& = 0), N = 20 x 20, 
run for 20000 steps (averages over 10 — 20 samples) The inset shows results for systems of size N = 10 x 10, 20 x 20, 40 x 40, 60 x 60 
at 5 — 0.1 from right to left. Right: Pf for different system sizes at fixed 8 = 0.1 at t — 5000. Dashed line in right panel is 
from numerical integration of the Master equation. 
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FIG. 9: (Colour on-line) Effects of interaction noise. S/N versus time t at different levels of external and interaction noise. 
Data is from simulations of systems with N — 10 x 10 agents. The left part shows 1? = (averages over 50 samples), the 
right one 1? = 2 (10 runs, data smoothened by running averages). F — 10, q — 100 in both panels. At large 7 no order is 
reached, independently of 5. For small enough 7 the system orders, the time-scale on which this occurs is largely reduced by 
the introduction of interaction noise, and increases as the interaction threshold is raised (note the different scaling on the time 
axes on the left and right). 



Oif — q,<7jf — 1. Distances different from zero or one thus indicate partial agreement between the two agents. 

In the following we will take the potential of two neighbouring agents to interact to be given by the following 
logit-rule |35j : 

- i + Jw (4) 

P is a control parameter allowing for the introduction of interaction noise. The case j3 = 00 here corresponds to the 
noise-free (zero temperature) case. If (3 — 00, agents with distance d > do are unable to interact, p(d > do) = 0, 
whereas interaction always occurs for pairs of agents with distance d < d a . d a is thus a threshold parameter, with 
large do corresponding to a regime of strong confidence of agents in other people's opinions, and small d to cases in 
which interaction is rare. In order to avoid confusion let us at this point stress that the role of the threshold do is 
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inverse to the one of i? in the nominal Axelrod model: large $ make interaction rare, whereas large do facilitate spin 
updates. 

Choosing finite values of (3 turns the hard threshold into a soft one. Interaction rates decrease smoothly with 
increasing distance. Crucially, at finite /?, interaction is always possible in principle, even if d > do- For (3 — finally, 
interaction is fully stochastic and independent of d. At any iteration, any chosen pair of neighbouring agents interacts 
with probability 1/2. 

Let us summarise the resulting dynamics: 

(i) Select one spin i S {1, ... , N} at random. Subsequently select one of its four nearest neighbours at random. 
Call this second spin j. 

(ii) Compute the Euclidean distance d(i,j) between i and j. 

(iii) If d(i,j) — both agents agree on all features. Interaction has no effect. If d = d(i,j) > then with probability 
p(d) as defined above spins i and j interact as in the nominal Axelrod model: one feature / is chosen at random 
so that a if ^ Oj/. Then set a if = Cjf- 

(iv) External noise. With probability 7 perform the following: choose one spin i and one feature / at random. Set 
cry to a value chosen randomly from {1, . . . , q}. 

(v) Resume at (i). 

A. Noise-free dynamics 

The behaviour of the noise-free system with Euclidean metric is described in Fig. [101 A transition between a 
disordered phase at low thresholds do and an ordered state at larger values of do is observed. The behaviour in these 
two phases is as follows: at low do only neighbouring agents with small differences in opinion can interact, so that 
the fraction of active bonds initially contained in the system is small. Dynamic arrest occurs quickly, and the system 
remains disordered. At large enough thresholds do the coarsening dynamics can persist until a fully ordered state is 
reached. Interestingly, as shown in Fig. [101 the critical value of the threshold do c does not depend much on the choice 
of F and q, and takes values do c ~ 0.4. The transition appears to become sharper at larger values of F (see Fig. ITU)) . 

Plotting Pp versus q at fixed F suggests that q plays no significant role in the Euclidean model. Only for small 
values of q can a dependence of Pp on q be detected. This invariance is intuitively to be expected as d is normalised to 
range between and 1 in the setup chosen here, so that q is merely a measure for how many discrete values can occur 
inbetween. Simulations with continuous opinions ranging in the interval [0, 1] (not shown here) reveal a behaviour 
very similar to the one depicted in Fig. 1101 We have also tested models with continuous opinions, in which both 
interacting agents agree on the average opinion of a given feature in case of interaction (with the same metric and 
kinetic constraints as before), and find similar behaviour as a function of do- Similar models are discussed in p^l[20|. 
mostly focussing on the case of one feature. 

B. Effects of external and interaction noise 

The behaviour of the Euclidean system under the influence of external and interaction noise is shown in Fig. 1111 
As seen in the main panel, a crossover between an ordered regime at low magnitudes 7 of the external noise and 
a disordered state at higher noise-amplitudes is found, very much like in the nominal Axelrod model. Interaction 
noise (finite (3) appears to have only little effect on this crossover for all values tested. Due to long equilibration 
times we have not performed a full analysis of the impact of external noise (7 > 0) in the large-A limit of the model 
zero-temperature ((3 — 00). The inset of Fig. [TTJ however demonstrates that at finite (3 the range of 7 in which the 
system orders is reduced as the system size is increased, similarly to what is found in the nominal Axelrod model. 
The ordering behaviour at small values of the external noise strength hence again appears to be present only in finite 
systems. 

C. Extremism and polarisation 

One particular potential application of metric Axelrod models is to study extremism, and the effects of a small 
number of highly polarised agents, e.g. with opinions at the extreme ends of the political spectrum. Having studied 
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FIG. 10: (Colour on-line) Pf versus interaction threshold do for the Euclidean model. Data is from simulations at different 
values of F and q (N = 20 x 20, averaged over 10 samples). Simulations are run until each bond has distance d = or d > do, 
i.e. until no active bonds are left in the system. 




FIG. 11: (Colour on-line) Effects of noise on the Euclidean model. Main panel shows Pf versus 7 from simulations of a 
system of size N = 20 x 20, run for 200000 macroscopic steps, averages over 10 samples (F = 10, q = 100). The interaction 
threshold is do = 0.3, i.e. the system is in the disordered phase if j3 = 00. Inset: Pf versus 7 at fixed J3 = 1 for system sizes 
TV = 40 2 , 20 2 , 10 2 from left to right (same simulation parameters as in main panel, 5 samples only for N = 40 2 ). 



the effects of uniform confidence thresholds on the opinion dynamics, we here extend the model to the case of 
heterogeneous, i.e. agent-specific thresholds. Here any agent i G {1, . . . , N} holds an individual threshold do^. In an 
interaction with a neighouring agent j he adjusts his opinion vector cr, only if d(i,j) < d,, i.e. do m Eq. (0]) is replaced 
by dcM- Agents with large do,i are thus likely to interact with others, and have a large tolerance against opinions in 
their surroundings. Agents who are unlikely to modify their own opinion vector are described by small interaction 
thresholds. Related work on other opinion dynamics models, mostly with continuous opinions and focusing on one 
feature or on nominal features, can be found in [HI, H2, HE HI] , see however also [2^. Our simulations here focus on 
the effects of extremists in the context of a multi-variate opinion dynamics model (F > 1) with discrete opinions and 
metric features. 

In this section we assume that the population of agents contains a fraction e of what we will call extremists. These 
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FIG. 12: (Colour on-line) Metric model with (left) and without extremists (right). In the left panel a fraction of s = 0.05 of 
extremists was used, with d e = 0.05. do = 0.6 as well as F — 10, q — 100 in both panels. Each panel shows the average opinion 
N^ 1 ^\ Oif on a given feature / as function of time for 10 individual runs of a system with N = 20 x 20 individuals, started 
from random initial conditions. Histograms of opinions at convergence are shown in the top part (data for 200 samples and 
from all F features used to generate histograms). . 

are agents whose opinion vectors take extreme values atf = 1 or aif = q at the beginning of the dynamics, and who 
are intolerant against other agents' opinion. In particular we assume an interaction threshold small d e for such agents. 
All other agents are taken to have a uniform threshold do as before and are initialised at random opinion vectors. The 
choices in the simulations presented in Fig. [12] are e — 0.05, d e = 0.05 and do = 0.6 (the latter threshold is chosen 
to ensure the system is in the 'active' ordered phase where dynamics persists long enough to prevent the system 
from remaining stuck in a configuration similar to the initial condition) . The figure shows the time evolution of the 
average opinion of the population on a given feature, along with the distribution of opinions at convergence. Apart 
from the initial conditions no stochasticity is present in the simulations shown in Fig. 1121 i.e. we have [3 — oo, 6 = 0. 
As shown in the left panel a small group of extremists can polarise the population, provided they are inert enough 
against adapting their own opinions. The figure shows 10 runs of the Axelrod dynamics. In each run extremists 
are chosen either to correspond to opinion states 1 or q. Extensions to two groups of extremists at either end of 
the political spectrum, or to cases in which extremists have polarised opinions only on some but not all features, 
are straightforward. The right panel shows the behaviour of the system in absence of extremists. Here the average 
opinion on each feature converges to a random value between 1 and q depending on the stochastic initial conditions, 
and the resulting histogram of final opinions is flat. 



VI. SUMMARY AND CONCLUSIONS 



We have extended Axelrod model for social influence to include varying interaction threshold, noise and metric 
features. In the basic model, which is two-dimensional, individuals are represented as multi-component spins and 
interact to become culturally closer, starting from a disordered random initial state. In our extensions some external 
and interaction noises are introduced, a confidence threshold limits the interaction, and a notion of distance between 
opinion is considered. We find that the confidence threshold does not influence the qualitative behaviour of the model, 
and that the typical transition, separating an ordered and a disordered equilibrium state, is preserved and determined 
basically by the initial probability of interaction between two individuals. The threshold limits the probability of 
interaction, and hence favors disorder. The introduction of an external noise brings in finite systems a continuous 
transition between an order-favoring and a disorder-favoring role of the noise, according to whether the noise is small 
or large respectively, independently of the threshold. This implies a discontinuity, in the region where the final state 
would be disordered in the absence of external noise. Such a discontinuity is removed in the thermodynamic limit. An 
interaction noise instead always favors order by reducing relaxation times, and does not alter the phase structure of 
the model. The other variant we study is the one with a notion of distance between opinions. We find that the model 
exhibits an order-disorder transition as the distance threshold is varied, consistently with the idea that the transition 
is the result of the competition between an ordering dynamics (the relative importance of which is determined by the 
distance threshold) and an initial disorder (as measured by size of the space from which the starting configuration 
is drawn at random). Moreover, the introduction of heterogeneous confidence thresholds in the context of metric 
Axelrod systems allows one to study e.g. the question of whether extremism can prevail in such models. We find 
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that the presence of a small fraction of individuals with a sufficiently rooted opinion can drive the whole population 
to the extreme ends of the opinion spectrum. Further application of heterogeneous interaction thresholds and metric 
features might include extensions addressing immigration or geographic barriers. Immigrants can for example be 
assumed to be more likely to interact with other immigrants than with members of the original population, and 
geographical barriers can be modelled by suppressing interaction at certain locations in space. This would lead to 
different interaction thresholds and tolerance levels, modulated either in space or dependent on the two agents picked 
for potential interaction. 
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